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We provide a detailed quantitative description of singular inflation. Its close analogy with finite¬ 
time future singularity which is associated to dark energy era is described. Calling and classifying the 
singularities of such inflation as finite-time cosmological singularities we investigate their occurrence, 
with special emphasis on the Type IV singularity. The study is performed in the context of a general 
non-canonical scalar-tensor theory. In addition, the impact of finite time singularities on the slow- 
roll parameters is also investigated. Particularly, we study three cases, in which the singularity 
occurs during the inflationary era, at the end, and also we study the case that the singularity 
occurs much more later than inflation ends. Using the obtained slow-roll parameters, for each case, 
we calculate explicitly the spectral index of primordial curvature perturbations n s , the associated 
running of the spectral index a 3 and of the scalar-to-tensor ratio r and compare the resulting values 
to the Planck and BICEP2 data. As we demonstrate, in some cases corresponding to the Type 
IV singularity, there might be the possibility of agreement with the observational data, when the 
singularity occurs at the end, or after inflation. However, absolute concordance of all observational 
indices is not achieved. On the contrary, if the singularity occurs during the inflationary era, this is 
catastrophic for the theory, since the observational indices become divergent. We also show how a 
Type IV singularity may be consistently accommodated in the Universe’s late time evolution, and 
we study the behavior of the effective equation of state. In addition, we investigate when inflation 
ends classically, in the context of our scalar-tensor model. Finally, we investigate which F(R) gravity 
can generate the Type IV singularity, with special emphasis on the behavior near the finite time 
singularity. 

PACS numbers: 04.50.Kd, 95.36.+X, 98.80.-k, 98.80.Cq,11.25.-w 


Introduction 

The appearance of singularities in classical cosmology is an inevitable inherent feature that almost all the theoretical 
descriptions carry along. In the Einstein-Hilbert description of local gravity effects, singularities are hidden behind the 
horizons of black holes, so the appearance of cosmological singularities at a large scale is rendered a rather unwanted 
feature of classical cosmological theories. Therefore, one of the most important quests in current cosmological research 
is to fully understand the presence of singularities, either these appear in the far past or in the far future. In the former 
case, the singularity is usually called the initial singularity, and a very fundamental question is posed, that is, whether 
the Universe really started in a singular state and continued its evolution entering an inflationary era |l| , or whether the 
Universe is described by some sort of a bounce [1J, in which case there is no initial singularity. The initial singularity is 
one point in spacetime at which spacetime is geodesically incomplete. This type of singularities is perfectly described 
by the Hawking-Penrose singularity theorems, given some time ago [3]. The geodesic incompleteness is practically 
ensured by the existence of time-like and null geodesics that cannot be extended to arbitrary values of their parameters. 
This constraint is equivalent to the general and stringent requirements of the Hawking-Penrose theorems, which are 
epitomized by the strong energy condition > 0, for every non-space-like vector K^, and also from the non 

existence of closed time-like curves. In addition, on these spacetime singularities, some of the curvature invariants 
diverge. However, sudden and Type II singularities that occur at a finite future time, can be quite milder from the 
initial spacetime singularity, in the sense that these are points at which it is possible that some observable quantities 
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(for example the Hubble rate or its higher derivatives) might become unbounded at some finite time, but this does 
not necessarily implies geodesic incompleteness. This type of singularities was studied by Barrow Q some time ago, 
and further developed in @41. According to Barrow’s description of the sudden singularities, at these singularities 
the strong energy and weak energy condition are satisfied, contraction is not required and the scale factor, along with 
its first derivative are finite. However, the second derivative of the scale factor and the pressure blow up at the finite 
time which these singularities correspond to. In addition, there exist sudden singularities of generalized type, at which 
the singular behavior occurs for higher derivatives of the scale factor. In principle, although that sudden singularities 
are scalar polynomial curvature singularities, these are much weaker than the Big Rip singularities, since geodesics 
are extendible through these points. In addition, these are not singularities of the crushing type, such as the Big Rip 
singularity I,®. 

Another consistent and important classification of future cosmological singularities, in which the value of observable 
quantities is taken into account, was given in [s], These singularities occur frequently in the context of modified gravity 
theories El- According to this classification, the observable quantities that are taken into account are the scale factor, 
the effective energy density, the effective pressure, the Hubble rate and its higher derivatives. It is obvious that some 
singularities of this type may not necessarily lead to geodesic incompleteness, thus can, in principle, be milder in some 
sense from the Universe’s initial singularity. The initial singularity and perhaps the Big Rip singularity are crushing 
singularities, in the following sense: at some time during the evolution, each point that lies on the spacelike three 
dimensional hypersurface t = const corresponds to a unique geodesic curve. Loosely speaking, these geodesic world 
lines form a congruence of geodesics that if these diverge in the past, or converge in the future to a unique point, then 
this point is essentially a singularity where the geodesic incompleteness occurs. 

Recently in Ref. @], Barrow and Graham studied a large family of finite-time cosmological singularities for a 
canonical scalar field with a power-law potential. Using a thorough qualitative analysis, they demonstrated that 
these weak finite-time singularities are an inherent feature of cosmological models with similar matter content. In 
addition, for appropriate initial conditions these models can describe large field inflation, with the difference that 
the physical system evolves to a singular state when inflation ends. Motivated by this work, we provide a detailed 
quantitative description of the singular inflation analysis, and classify the singularities found in [§}, according to the 
cosmological singularities classification given in 0-0 SI • As we shall demonstrate, there is some difference in our 
analysis, in comparison to the one presented in J8| , traced in the case at which higher derivatives of the Hubble rate 
are taken into account. This feature characterizes the Type IV singularities, according to the finite time cosmological 
singularities classification given in @j. More importantly, as we explicitly show, the existence of a Type IV singularity 
in an inflationary model can have dramatic effects on the slow-roll inflation parameters, rendering them singular at 
finite time, if the singularity occurs during inflation. We have to note that it is possible the slow-roll inflation does 
not work near the Type IV singularity. Indeed, the corresponding approximation may be broken especially, if the 
singularity occurs during inflation. Moreover, if inflation ends exactly at a Type IV singularity, this can be compatible 
with current observation data, and interestingly enough, after the Type IV singularity, the Universe is accelerated in 
a quintessential way. In addition, we study the evolution of the effective equation of state as a function of cosmic 
time and we also address the question when inflation ends, providing a detailed quantitative analysis. With regards 
to the question when inflation ends, we found a particularly interesting solution which occurs in the case of a Type 
IV singularity, according to the finite time cosmological singularities classification provided in @]. As we show, there 
appears the physically appealing picture that, if certain requirements are met, inflation ends, classically at least, at 
a non-singular state at a time t e , before the Type IV singularity occurs, and the Type IV singularity is reached in 
a future time, t s , with t s t e . Thus in a Universe that evolves in time with a Hubble rate that has a Type IV 
singularity, inflation can end in a non-singular state. 

In the study of future Universe after ACDM-like epoch, the finite-time future singularities that may occur, are 
classified according to Ref. @|. However, most attention was given to Type I (Big Rip) future singularity (for 
incomplete list of references on that, see Refs. jTol. fl2|) or Type II singularity, commonly known as a sudden future 
singularity (for an incomplete list of references on that, see Ref. 0)- These future singularities may occur for phantom, 
or quintessence-like, dark energy Universe, even if currently it shows nearly ACDM-like evolution. As the classification 
of these singularities is applied to accelerating dark energy Universe, it is evident that it may be applied in the same 
sense to the classification of singularities after singular inflation. Hence, in this paper we adopt just the same notations 
for all four types of singularities, even after inflation. Of course, one should bear in mind that it is not physically the 
same as future singularities after ACDM epoch, because finite-time singularity may occur prior to dark energy epoch 
or after current dark energy epoch. In addition, it may also occur after dark energy epoch. In this respect, we shall 
focus on the Type IV future singularity, because it is expected that Universe may survive after passing through such 
singularity smoothly, since it is not a crushing type of singularity. Thereby, it could be possible that between (singular 
perhaps) inflation and the ACDM epoch, the Universe passed through type IV singularity. Furthermore, for unified 
inflation-dark energy models of the sort introduced in Ref. [13] , we may get the following really physical appealing 
possibility. The future singularity of Type II, III or IV, may occur just because of the epoch of singular inflation, so 
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that dark energy era could be de Sitter-like or pure ACDM. In other words, the future singularity may be predicted 
by only the study of unified inflation-dark energy cosmology, like in the examples of Ref. [13l |. However, by studying 
only the dark energy epoch, one discovers that it is ACDM-like and therefore, the future universe is regular. Thereby, 
based only on the analysis of dark energy epoch, one may arrive to wrong conclusions about the future universe. We 
shall exemplify a scenario like this in a section, later on in this paper. 

Finally, owing to the fact that the finite time cosmological singularities classified according to Refs. Q, frequently 
occur in the context of modified gravity, we investigate which F(R) gravity can generate the Type IV singularity, 
providing the exact form of the F(R) gravity near the Type IV singularity. As is expected, the resulting F(R) gravity 
is approximately of polynomial type near the singularity, which can be fine-tuned to be of the Einstein-Hilbert form 
gravity plus curvature corrections. 

This paper is organized as follows: In section I, we quantify Barrow’s singular inflation description, in the context 
of scalar-tensor gravity, by using very well known scalar-tensor reconstruction techniques EMI. Using very general 
arguments, we thoroughly describe all the singularity types that can occur in these theories, and we exemplify our 
results using illustrative examples. In addition, we investigate the effects of a singularity in the slow-roll inflation 
parameters and also study the behavior of the effective equation of state parameter as a function of the cosmic time. 
Moreover, the question when inflation classically ends is addressed too. In section II, we investigate a physically 
appealing scenario, according to which, in the end of inflation, the Universe experiences a Type IV singularity and in 
addition the dark energy era is quintessential or approximately de Sitter and at finite time after the dark energy era, 
the Universe experiences again a Type IV singularity. In section III we investigate which F(R) gravity can realize the 
Type IV singularity, which proves to be the most relevant type of singularities. We address this issue using qualitative 
arguments and we quantify our study by using very well known reconstruction techniques, in order to find the F(R) 
gravity near the Type IV singularity. In section IV we briefly investigate the impact of the finite time singularities 
on the long-wavelength scalar perturbation modes, with special emphasis in the Type IV singularity. In section V 
we discuss the graceful exit problem and this is affected by the presence of finite time singularities. Our conclusions 
along with a brief critical discussion, follow in the end of paper. 


I. SCALAR-TENSOR GRAVITY ANALYSIS 


In this section we shall provide a quantitative analysis of the singular inflation physics and its consequences on 
observable data, in the context of a general scalar-tensor gravity theory containing one scalar field. For related studies, 
the reader is referred to Refs. I9L [Taj. 

We start off our analysis by recalling how the finite time cosmological singularities are classified according to 
Refs. Q. For a detailed presentation of these issues consult also jlfil j. As was shown in Ref. sm , the finite-time 
future singularities are classified as follows, 

• Type I (“Big Rip”) : When t -A t s , the scale factor diverges a, the effective energy density p e g, the effective 
pressure p e g diverge, a —> oo, p e g —> oo, and beff| —> oo. This type of singularity was presented in Ref. 0 and 
studied in Ref. Q. 


Type II (“sudden”) EJ: When t —» t s , the scale factor and the effective energy density is finite, a 
but the effective pressure diverges \p e g\ —> oo. 


a s , Peg ^ Ps 


• Type III : When t — 
pressure diverge, p e g 


t s , the scale factor is finite, 


a s but the effective energy density and the effective 


OO, beff| -A OO. 


• Type IV : For t —> t s , the scale factor, the effective energy density, and the effective pressure are finite, that is, 
a -A a s , Peff A► p s , beffl -A Ps , but the higher derivatives of the Hubble rate H = a/a diverge. 

The Type III and Type IV singularities were extensively studied in Ref. 9j. Notice that, the effective energy density 
p e s and the effective pressure p e g are defined in the following way, 


Peg = —gH 2, 


Peg = 


b( 2ff+3H b 


(1) 


In Ref. (Hi, a scalar field model realizing inflation with a Type IV singularity has been proposed. In this section, by 
using the general formulation of the scalar reconstruction of Refs. 13j, 141, we construct general scalar-tensor models 
that realize Type II and IV singularities. We consider the following non-canonical scalar-tensor action, which describes 
a general single scalar field model, 


S = 


d A Xy/~Fg | — -V(<t>) + Tmatterj • 


(2) 
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The function ui((f>) is called the kinetic function, V{(j>) is the scalar potential, and both are assumed to be functions 
of the scalar field </>. In addition, we assume a flat Friedmann Robertson Walker (FRW) metric of the form, 

ds 2 = — df 2 + a 2 (t) d.T 2 . (3) 

i 

The function w((/>) in Eq. © is irrelevant and we can therefore absorb this function by redefining the scalar field <f>, 
in the following way, 

<4> = J d<j>y/u(<t>) , (4) 

where we assumed that w(</>) > 0. Then the kinetic term of the scalar field in the action ([2]) is rewritten as, 

- u{<j>)d tl <j>d li <l> = . (5) 

It is convenient for later purposes, if we keep w(</>), although we can absorb into the redefinition of the scalar 
field. In the case that < 0, the scalar field becomes ghost corresponding to phantom dark energy and in this 
case, instead of (01, by using the redefinition, 


we find that, 


<P = 



d(j)y/—ui(<t>) , 


( 6 ) 


- i9 m v ?( 9 / V, (7) 

instead of the expression given in Eq. © . In the case that the scalar field is a ghost field, the energy density becomes 
unbounded from below in the classical theory, but in the quantum theory, the energy becomes always bounded from 
below, with the cost that there appears a negative norm. 

From the action of Eq. ©, we obtain the following expressions of the energy density and the pressure : 

P = ^w(</>)^ 2 + V((/>), p= ^u;(»</> 2 - V{<j >). (8) 

Therefore, the scalar potential V(</>) and the parameter w(</>) can be written in terms of the Hubble rate and its 
derivatives as follows, 


c (</»)</> 2 = ~ H , V (</>) = 4 ( sH 2 + H ) ■ ( 9 ) 

Consider the following solution, 

4> = t, H = (10) 

If we assume that w{(f>) and are given by using the single function /(</>), as follows, 

aw = ~m, v{<t>) = 4 (3/(<^) 2 + rw) , (ii) 

we find that the solution ca is indeed a solution of the FRW equations of the system. 

It can be easily verified that the equation which results by the variation of the action, with respect to </>, namely, 

0 = w{<f>)<j> + {4>)(j) 2 + 3 Hu}(4>)4> + V'((j>) , (12) 

is also satisfied by the solution (fTOl) . In effect, the arbitrary Universe’s evolution, expressed by H(t) = f(t), can be 
realized by an appropriate choice of ui(<j)) and V(</>). In other words, by defining the particular type of the Universe’s 
evolution, the corresponding scalar-Einstein gravity may be found. 

In the case of the Type II and IV singularities, the Hubble rate H(t) may be chosen in the following form: 


H{t) = hit) + f 2 (t) [t s - t) a . 


(13) 
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Here fi{t) and f 2 (t) are smooth (differentiable) functions of t and a is a constant. If 0 < a < 1, there appears Type 
II singularity and if a is larger than 1 and not integer, there appears Type IV singularity. When a is given by two 
integers n and m by 


n 

2m + 1 


(14) 


we can consider the region where t > t s . 

Then, by using CD. we find 

w(</>) = “ ^ {/lO) + /2WO % - </>)“ - af 2 ((j>) (t s - </>) Q_1 | , 

v(<t>) = 13 (/ i (^) + f 2 (4>) ( t s - 4>T) 2 + /((</>) + f' 2 (4>) (ts - <t>) a - oif 2 ((j>) (t s - </>) a_1 |, ( 15 ) 

which determine a general model which generates the Type IV singularity when a > 1. When t = <j> ~ t s , if a > 1, 

the expression of w(<j>) in (fl5l) gives 

u(4>) - \f[(t s ). ( 16 ) 

Hence, although the scalar field is canonical if f[(t s ) < 0, the scalar field becomes ghost, corresponding to the phantom 
dark energy if f[(t s ) > 0. In other words, the Type IV singularity is realized not only in the case that the scalar 
field is canonical, but also the singularity can be realized even in the phantom phase, where the phantom dark energy 
dominates. 

We first consider the simple case that f\(t) = 0 and f 2 (t) = fo with a positive constant /o- In the neighborhood of 
t = t s , we find that, 

w(0) 0) 0 " 1 , VW ~ (t s - ^)“" 1 , (17) 

and by using d4j), we find 


272^/0 

k (a + 1) 


(t s 



Consequently, the scalar potential reads, 


V{p) 


otfo 

K 2 


«(a + l)^^ “+ 

2\/2a7o 


(18) 


(19) 


Therefore, when the following condition holds true, 


-2 < 


2(a~l) 

a + 1 


< 0 , 


( 20 ) 


there occurs the Type II singularity. Accordingly, the Type IV singularity occurs when the following holds true, 


0 < 


2(a-l) 

a + 1 


< 2 . 


The above results seem to be a little bit different from those appearing in 
potential is of the form V(<p) ~ ip n an d therefore we may identify, 


( 21 ) 

As was qualitatively shown in Q, the 


n = 


2(q-i) 

a + 1 


( 22 ) 


In the claim of Q, when 1 > n > 0, the Type IV singularity could occur, in which case H diverges, and when 
k + 1 > n > k, the fc-th derivative of H, that is H^ k+2 \ diverges. Then the Type IV singularity should occur, even 
if n > 2, which is in conflict with Eq. (ED- In the qualitative analysis of Q, it has been assumed that (p ^ 0 at 
t = t s when 3 > a > 2, that is, H is finite, but H^ diverges at t = t s , but as we can find from (|18[) . ip behaves as 

ip ~ (t s — t) 2 and therefore </> vanishes at t = t s . This shows that the arguments in |8| might not apply in our case, 
and specifically when 3 > a > 2 or equivalently n > 2. In addition, the Type IV singularity does not occur when 
n > 2. 

We should note that in Eq. urn 
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• a < — 1 corresponds to the Type I singularity. 

• — 1 < a < 0 corresponds to Type III singularity. 

• 0 < a < 1 corresponds to Type II singularity. 

• a > 1 corresponds to Type IV singularity. 

In the case that A (A = 0 and A (A = /o, we have assumed that fo should be a positive constant in order for the 
Hubble rate H given in Eq. (fl3l) to be positive. Then in the case that a is negative, which corresponds to the Type I 
and II singularities, w(0) in Eq. (1171) and therefore the scalar field becomes ghost corresponding to a phantom era. In 
the case that a is positive, which corresponds to the Type II and IV singularities, the potential V(ip) in (fl7l) becomes 
unbounded from below, but the expression in (1171) is only valid at t ~ t s and the full expression in (1151) gives, 

V(<t>) = {3/o (t s - tf a - af 0 (t s - tf- 1 } , (23) 

which is bounded from below if a > 1. Even if a < 1, by choosing fi(cj>) and / 2 (</>) properly, we can make V(ip) 
bounded from below and even to be positive semi-definite. In the case that fi(t) = 0 and f 2 (t) = fo , the expression 
n = implies that, 

• n > 2 corresponds to the Type I singularity. 

• —2 < n < 1 corresponds to Type III singularity. 

• n < —2 corresponds to Type II singularity. 


• 1 < n < 2 corresponds to Type IV singularity. 

Even in case that fi(t) = 0 and f 2 (t) = fo, with fo a positive constant, in Eq. (fI71) . the slow-roll inflation could 
occur, and at the end of the inflationary era there could occur the Type IV singularity. We shall elaborate on this 
further later on in this section. 

We may consider more general models, by choosing a more general form for the arbitrary functions fi(t) and / 2 (A- 
In this way, we can realize almost all kind of the cosmologies possessing a Type IV singularity at t = t s . For example, 
if we choose, 


or more generally, 


fi(t) 


A 

\Ao + A 


A (A 


fit 2 

4 + 1 4 ’ 


A (A 


A 

Ao+*")" ’ 


A (A 


ht 2 

4 + t 4 ’ 


(24) 


(25) 


the Universe becomes a de Sitter spacetime, with the Hubble rate H becoming a constant, that is, H ~ in the 
early Universe. In addition, H behaves as H ~ jA at late time, although the Type IV singularity appears at t = t s . 
If we choose /2 small enough, the second term (1131) can always be neglected, only if we consider the development in 
the expansion of the Universe. This argument is quite general, because we can always choose the second term to be 
small enough, compared to the first term in (fl3l) . Therefore, the Type IV singularity is not always relevant to the 
expansion history of the Universe. However, as we will now explicitly demonstrate, the Type IV singularity affects 
strongly the slow-roll parameters, if the singularity occurs during the inflationary era. 

The slow-roll parameters e, 77 and £ by H can be expressed as follows EHOSEI, 


e = 


77 = 


f = 


= _L (fy \ 2 (VW \ 2 = 2 _i_ (V'( 4 >)\ 


2k, 2 \dip 

1 

k 2 U(<A 

V'tt) 

k a V ( 4 >) 2 u }((/)) 2 


V u(</>) J 2k 2 w (0) V v(4>) J 


^4 ( 4 ) v W + (it 


dip dtp \dipj 

u"(4>) 
2 u)(<t>) 


u'(4>) 

w(<A 


v"(4>) 


y'(4>) - 


K 2 V (( j )) 


u'(4>) 

2 ui (( j )) 2 


V(4>) 


w(</>) 


v"(4>) 


3 ijj'( 4>) 
2 u>(4>) 


V'W + V'tt) 


(26) 


We should note that the derivatives of V(<j>) and u>(<j>) diverge at (f = t = t s in general. 
There could be three cases, 
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1. The Type IV singularity occurs during the inflationary era. 

2. The inflationary era ends with the Type IV singularity. 

3. The Type IV singularity occurs after the inflationary era. 

We first consider the case [T] that is, the case in which the singularity occurs during inflation. Then, the above 
slow-roll parameters may diverge, which may be in conflict with the current observational data. In the same line of 
research, there are several cases, in which the aforementioned argument holds true, for example, 

• When fi{t s ) = f[(t s ) = 0, V(<j>) and behave as (t s — f )“ -1 near the singularity as in (flTl) . Therefore we 


find, e, r] ~ (t s - t)~ 


e ~ Us - 1) 


—a—3 


and therefore all the slow-roll parameters diverge at t = t s 


• When fiit s ) / 0 but f[(t s ) = 0, <*;(</>) behaves as (t s — t) a , again but V{<f>) is finite if a > 1. Furthermore if 
2 > a > 1, V'{4>) ~ its — t) a ~ 2 and therefore e ~ (f s — t ) a ~ 3 , ~ (t s — 1)~ 2 , £ 2 ~ (t s — t)~ 4 , and therefore all 

the slow-roll parameters diverge at t = t s . 

All the parameters e, rj, and £ 2 are finite for a > 4, and also when, fi{t s ), f[(t s ), fiit s ), f[ 3 \t s ), and f[ 4 \t s ) do 
not vanish. In this case, however, these parameters do not depend on f- 2 it s ) and its derivatives, and consequently the 
slow-roll parameters, do not depend on the singularity. 

In all the above considerations, we have assumed that the Type IV singularity may occur during the inflationary 
era, as in case CD We may, of course, as in case [2] or [3] consider models that the singularity could occur after the 
inflationary era. If there is any matter coupled with the higher derivative of the curvature, which includes the higher 
derivatives of H, the Type IV singularity may affect the matter content, but owing to the fact that such a coupling 
is not renormalizable, this could be relevant only at a very high energy region. 

Let us now address in detail the remaining cases [2] or [3] We start off with case[ 2 l in which case inflation is assumed 
to end at the time tf which is equal to tf = t s , with t s . which appears in the Hubble rate of Eq. (1131) . being the 
time at which the possible singularity may occur. We also assume that /i(t) = 0 and / 2 (f) = /o, so that the Hubble 
parameter and the corresponding scalar potential near the singularity t ~ t s , are equal to, 


Hit) = f 0 its - t) a , V{4>) = -af 0 its - </>) 


a-1 


(27) 


It is worth working using the canonical scalar field <p, instead of <j>, in order to make direct contact with the model of 
Ref. [§]. In terms of the canonical field <p, the scalar potential near the singularity is equal to, 


with n given in Eq. 


In addition, Vo in Eq. (1281) . is equal to 


V(<p) * V 0 <f n 

, but we quote it again for simplicity, 

a + 1 


Vn = — 


afo 


2(a — 1) 

re (a + 1 ) V “ +1 


2v^/o 

In terms of n, the Hubble rate is written as follows, 

Hit) = f 0 its - t)~^+ 


(28) 


(29) 


(30) 


(31) 


We shall investigate the impact of the fact that tf = t s . on the inflation slow-roll parameters and consequently to 
the observational indices corresponding to inflation. The slow-roll parameters we shall study are the ones given in 
Eq. (1261) . which when defined in terms of the canonical scalar field tp are equal to, 


= _L (Y1M) 

2re 2 'v Vip) ) 


1 

V = — 

if * 1 


V"(<p) 

Vip) 


e = -o 


1 (V'{p)V"'ip) 


V{pf 


(32) 


with the prime this time denoting differentiation with respect to p. It is worth to recall at this point how the slow- 
roll parameters are derived for the case of a canonical scalar field. For a detailed account on these issues, consult 
Refs. [il. Il7j|. In general, the slow-roll condition relies upon the constraint, 


1 


< V ip ). 


(33) 
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and that this constraint holds true for an extended period of time. This condition ensures a long and finite acceleration 
era and it is known as the first slow-roll condition. In order for the first slow-roll condition to last over an extended 
period of time, a second condition must imposed, 


\(p\ < 


dv{y) 

dp 


(34) 


which is known as the second slow-roll condition. In virtue of the equation of motion satisfied by the canonical scalar 
field in a flat FRW background, the condition (1511) can be written as, 


1 $ < ZH\(p\. 


(35) 


With the two slow-roll conditions holding true, the canonical scalar field equations for the scalar field reads, 


1 dV(p) 
^ 3 H dp 


while the FRW equations become, 


3 H 2 ~ 8 t jGV(p). 


(36) 


(37) 


It can be shown that the first and second slow-roll conditions, namely Eqs. (1511) and (1351) . are equivalent to the 
following relations, 


V'{<P) 

V(<p) 


< 2k 2 , 



(38) 


with the prime denoting differentiation with respect to the canonical scalar field. The above equations can be written 
as follows, 


e « 1 , r) < 1 , 


where e and 77 are the slow-roll parameters, 

= J_ (YM) = 1 

2ft 2 \ V(p) J ’ ft 2 VR((/9 2 ) ) 


(39) 


(40) 


Notice that the conditions (1551) are actually the slow-roll conditions. 

The observational quantities we shall study are the spectral index of the primordial curvature fluctuations n s , the 
tensor-to-scalar ratio r, the associated running of the spectral index a s , which are defined in terms of the slow-roll 
parameters as follows [lj, 1131 j 


n s ~ 1 — 6e + 2rj , r = 16e , a s ~ 16e?7 — 24e 2 — 2^ 2 . (41) 

The latest Planck observational data [l8| constrain the above observational quantities as follows, 

n s = 0.9603 ±0.0073, r < 0.11, a s = -0.0134 ± 0.009, (42) 

and the BICEP2 [li| observations strongly indicate that, 

r = 0.2t° o ° o l (43) 

Note that the Planck data are currently considered to be more reliable than the BICEP2 results. Finally, the values 
of the scalar field at the beginning and at the end of inflation, and pf correspondingly, are chosen is such a way 
so that the Lyth bound [2Cj is respected, so that, 


\npi - np f \ > 1, (44) 

when the BICEP2 results are taken into account, while consistency with the Planck data requires, 


\npi - Kpf\ < 1 . 


(45) 
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Practically, these requirements imply that the number of e-folds until inflation ends, is N ~ 50 for the BICEP2 result, 
while for the Planck, N ~ 60. Having all the above at hand, we can perform a full quantitative analysis of the 
observational implications that a Hubble rate and the corresponding scalar potential of the form given in Eq. m 
would have. A similar detailed analysis for the case of chaotic inflation potentials was performed in [23|. For the 
potential of Eq. m, the slow-roll parameters are equal to, 


e(if) = 


2k 2 (t s - p) 2 1 


r](ip) = 


(—1 + n)n 
« 2 (U - p) 2 ’ 


2 _ (-2 + n)(—1 + n)n 2 
K 4 (t s - p) A 


(46) 


So let us now investigate what will happen if tf = t s . Under the slow-roll approximation, inflation will actually 
end, when the scalar field acquires it’s final value < pf 1 at which the slow-roll parameter e(pf) (or r](p))) will be of the 
order e(pf) ~ 1. From Eq. (H51) . this would imply for tpf that, 


t s -tpf = 


n 

\[2k 


(47) 


In order to have a clear picture for the ending of inflation, we shall express the slow-roll parameters as a function of 
the number of e-folds until the end of inflation, which is defined to be equal to, 


N = — 


H(t)dt , 


(48) 


and for the Hubble rate given in Eq. m, this is equal to, 

N = —(—2 + n)(—t + t s ) _2 + n — —(—2 + n)(—tf + t s ) _2 + n , (49) 

which is clearly divergent for tf = t s , if n > 2. Recall that the e-folding number is constructed in such a way so that, 
when inflation ends, it has to be exactly equal to zero, with it’s value increasing as p tends to it’s initial value pi. 
Therefore, we find a first point of inconsistency caused for some values of n. Recall from our previous classification 
of singularities that the n > 2 case corresponds to the Big Rip singularity, so in the case of a Big Rip singularity, we 
come across to a severe divergence of the e-folding number. So the Big Rip case is an idiomorphic case to be excluded 
from our study, at least when case [2] is considered. Therefore, in the following we shall assume that n < 2, at least in 
the context of case[2j 

In the slow-roll approximation, the e-folding number N, may be approximated by the following relation 


N ~ k 2 




V(<P) 

V'(p) 


d p, 


which for the scalar potential of Eq. (E51) . reads, 


N ~ k 2 - 


tsP P_ t s pf 

n 2 n n 


(50) 


(51) 


Substituting the value of pf from Eq. (1471) in the above equation, the e-folding number becomes, 

n 2 — 2 n 2 (t s — p) 2 


N~ — 


An 


so by solving the above equation with respect to t s — p we obtain, 


t s -p = 


ANn - 


2k 2 


(52) 


(53) 


We can substitute t s — p appearing in Eq. (l53l) . to the slow-roll parameters of Eq. (l46l) . and finally we will have the 
slow-roll parameters as functions of the e-folding number for the whole inflationary era. By doing that we obtain, 


4iV + n ’ 


2(—1 + n) ^2 $(—2 + 7l)(—1 + 77.) 


AN + n ’ 


(AN + n) 2 


(54) 
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Observational Indices 

Spectral Index n s 

Tensor-to-Scalar ratio r 

Running of Spectral Index a s 

TV = 50, n = 1.9 

n s ~ 0.961367 

r ~ 0.15057 

a s ~ -0.000765389 

TV = 50, n = 1.382 

n s — 0.966412 

r ~ 0.109801 

a s ~ -0.000727579 

TV = 50, n = 0.5 

n s ~ 0.975062 

r ~ 0.0399002 

a s ~ -0.000497509 

TV = 54, n = 1.48 

n s ~ 0.967997 

r ~ 0.108884 

a s ~ -0.000588614 

TV = 54, n = 1.8 

n s — 0.965106 

r ~ 0.132231 

a s ~ -0.000640852 

TV = 54, n = 0.9 

n s ~ 0.97326 

r ~ 0.06639 

a s ~ -0.000493138 

> 

cn 

o 

3 

II 

i— 1 

O 

n s ~ 0.969384 

r ~ 0.112536 

a s ~ -0.000506685 

TV = 60, n = 1.1 

n s ~ 0.974285 

r ~ 0.0729988 

a s ~ -0.000426636 

TV = 60, n — 0.9 

n s ~ 0.978387 

r ~ 0.0399002 

a s ~ -0.000359312 


TABLE I: Values of the spectral index of primordial curvature perturbations n s , the associated running of the spectral index 
a B and of the scalar-to-tensor ratio r, for various values of the e-folding number TV and of the parameter n, for the case that 
the singularity occurs exactly at the end of inflation (Case [2]). As it can be seen, the Type IV singularity gives results that are 
compatible to BICEP2 and Planck data, when the spectral index n s and the scalar-to-tensor ratio r is considered 


so that finally, the observational indices read, 

4(7V — 1) — n 

ria — - 

47V + n 


r — 


16n 

47V + n 


a s — — 


8(n + 2) 
(47V + n) 2 


(55) 


Having these at hand, we can promptly examine which values of n, are ruled out by observations, by taking the e- 
folding number to be TV = 60, TV = 50 or TV = 54. By looking at the Tabled] it is obvious that, for (TV, n) = (50,1.382), 
(TV, n) = (54,1.48) and (TV, n) = (54,1.8), the results are compatible with the BICEP2 observational data, with regards 
to the spectral index n s and the scalar-to-tensor ratio r. However, in all the cases, the running of the spectral index 
is smaller than the Planck one. For (TV, n) = (50,1.382) and (TV, n) = (54,1.48), the spectral index n s and the ratio 
r are compatible to the Planck observational data. Recall that when 1 < n < 2, the physical system develops a Type 
IV singularity and by looking Tabled] it is obvious that, if inflation ends at a Type IV singularity, the results can be 
compatible to the Planck and BICEP2 observational data. 

Now we consider case |3] in which case t s is considered to occur much more later than the inflation ending time 
tf, that is, t s tf. Then, contrary to case [2] there is no inconsistency in the definition of the e-folding number in 
Eq. (1491) . since it is always finite, and in addition it is equal to zero when inflation ends, as it should be. So practically 
all the values of n are allowed and all singularities may occur at late times, with no inconsistency. Let us see how 
the observational indices change in this case. Following the lines of research of our previous analysis, in the slow-roll 
approximation, inflation ends at, e(tpf) ~ 1, so in this case we have, 




(~ l + n) 
\/2k. 


so in the slow-roll approximation, by substituting iff from Eq. (1561) . the e-folding number TV reads, 


TV- - 


—2\/2 K(tf — t s )n + n 2 + 2« 2 (f/ — <p){tf — 2 t s + tp) 


4 n 


and by solving with respect to p, we get, 


<P = 


2 ^4ft 2 f s + ^ 1 16 kH 2 — 8k 2 (^ — 1 — 47Vn — 2 n 2 t 2 + 4 n 2 tft s + 2a + 2\Z2n(tf — t 8 )n — (n + l) 2 ^ ■ 


(56) 


(57) 


(58) 


Substituting the above expression in Eq. (1461) . we obtain the slow-roll parameters, with e being equal to, 


1 + 47Vn + (n + l) 2 — 2 (n + 1) (l + \J2 (tf — t s )n ) + 2 (tf — t s )n (>/2 + tfn — t s K) 


( 59 ) 
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Accordingly, 77 is equal to, 


^ (l + 4iVn + (n — l ) 1 2 — 2 (n 
while, £ 2 reads, 


f 2 


(l + 4An + (n — l ) 2 — 2(n 


2 (n — 1 )n 

1) (l + V2(t f — t s )«) + 2 (tf — t s )n (V2 + tfK — t s n )) 

4 (n — 2)(n — l)n 2 

1) (l + \/2(t/ — t s )«) + 2 (tf — f s )ft (>/2 + tfK - t s n )) 


Having the slow-roll indices at hand, we can promptly calculate the observational indices, which are, 

4 (n — 1 )n + (41V n — 5?r 2 — 2\/2 (tf — t s )nn + 2 (tf — t s ) 2 n 2 ) 

j) ~_____ 

(l + ANn + (n + l) 2 — 2(n + 1) (l + y/2(tf — t s )«) + 2(tf — t s )n (a/2 + tfK — Uk)) 

16n 2 

1 + 41Vn + (n + l) 2 — 2(n + 1) (l + \[2 (tf — t s )tij + 2 (tf — t s )n (-\/2 + tfK — t s n ) 

8n 2 

d s — y 

(l + 41Vn + (n + l) 2 — 2(n + 1) (l + y/2(tf - t s )«;) + 2(1/ — t s )n (-\/2 + tfK — 1 s k)) 
x ^6 — 5(n + 1) + (n + l) 2 + ^1 — 8 Nn + (n + l) 2 + (n + 1) ^—2 + 4\/2(1/ — 1 s )k^ 

—4(1/ — 1 s )k ^-\/2 + tfK — t s iij ^ . 


(60) 


(61) 


(62) 


The above expressions are quite complicated, but can be simplified if certain assumptions hold true. It is quite natural 
to assume that t s is a quite large number, so that the singularity occurs much more later than the present time. But 
even if we choose l s to be equal to the present age of the Universe, t s is quite large, so let t s = 6.636 x 10 51 GeV _1 , 
which is the present age of the Universe. In addition, k = 8ttG = 2.0944 x 10 - 18 GeV -1 , with G, Newton’s gravitational 
constant. So obviously, the term 2 k 2 (1/ — t s ) 2 dominates (it is of the order 10 68 ) in all relations of Eq. (1621) . so the 
observational indices may be simplified as follows, 


n s ~ 1, 


256 


2 K 2 (tf - t s ) 2 
8 n 2 


(fc 2 (l/-l s ) 2 ) 


< 1 , 

< 1 , 


(63) 


with a s > 0. Therefore, the Planck constraints are not satisfied, when the spectral index n s and the running of the 
spectral index a s are taken into account. Now a question rises, what if the time that the singularity occurs, occurred 
in the past, but after inflation? In that case, it is possible that the observational constraints are satisfied, for example, 
if t s ~ 2 tf, that is, the singularity occurs after the end of inflation, for ( N,n ) = (55,1.55), we get, n s ~ 0.967953 and 
r ~ 0.108352, which are compatible with Planck and BICEP2 results. 

In conclusion, we can say that the Type IV singularity offers a good phenomenology, since in all the cases that 
n contributes in the resulting expressions of the observational indices, some of the values with 1 < n < 2 , result to 
observational indices compatible with observations. 

Coming back to the paradigms for the functions /i(t) and / 2 (f), instead of choosing (1M)) . we may choose, 


h{t) 


h 

\Ao + 


+ fo ^ 


f2(t) 


f2t 2 

4+t 4 ’ 


(64) 


where it is assumed that ^ f 0 > 0 and also that / 2 is small enough. Then, the behavior of the early time is not 

different from that given in Eq. (l24l) and we find H ~ ^ but even at late time (t —> 00 ), the Hubble rate becomes 
a constant H ~ fo, which may describe the dark energy in the present Universe. Then the model unifies inflation in 
the early Universe with the accelerating expansion in the present Universe. If we choose to large enough, the Type 
IV singularity may occur in the Universe where the dark energy dominates. 

We now study the effective equation of state (EoS) parameter w e s, which is equal to, 


r . 

W e ft = - = -1 ~ 

P 


2 H 
3H 2 ' 


( 65 ) 
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This parameter will reveal whether the scalar field will go to phantom or not, during the process of cosmological 
evolution. We shall be mainly interested for the model with Hubble rate given in Eq. m, with fi{t) = 0 and 
/ 2 (f) = /o > 0. For this Hubble parameter, the EoS parameter reads, 


WeS = -1 - 


2a(—t + t s ) 1 “ 

W 


( 66 ) 


Obviously, the EoS parameter can be singular if a > — 1, and in the converse case it is non-singular. Let us investigate 
if the EoS evolves to phantom or not. When a < — 1 and since fo > 0, the sign of the fraction appearing in Eq. (EH) 
is determined by the actual value of —a — 1, for t > t s . We shall consider cosmological times after the singularity, 
that is, t > t s and in addition, a is assumed to be of the form given in Eq. m, with n an even integer. Then, when 
t > t s , and regardless if a > — 1 or a < — 1, the EoS describes an quintessence accelerating Universe, since, the EoS is 
always w e s > —1. In the case that a < —1, the EoS is regular at t = t s , while for a > —1, it is always singular. Let 
us make contact with the singularity study we performed earlier. When a > —1, the EoS is always singular at the 
transition point t = t s . For a cosmological time with t > t s , then the EoS is u> e ff > —1, so non-phantom, quintessence 
acceleration occurs. Notice that when a > 1, the singularity is Type IV, so the fact that the EoS is always w e g > — 1, 
has very appealing physical consequences, as we will demonstrate shortly. 

In conclusion, the EoS of state analysis revealed a physically appealing case which occurs when t s is considered to 
be an early time in the evolution of the Universe. Consequently, a natural question to ask is, is it possible for inflation 
to end at this Type IV time and more generally, when does inflation ends. As we explicitly demonstrated earlier, 
if inflation ends exactly at a Type IV singularity, which occurs at t = t s , it may be possible for the theory to be 
compatible with the observations. Having these issues in mind, we now perform a qualitative analysis of the inflation 
ending time. A similar analysis on the inflation ending was performed in jlSl j. Practically, inflation ends when a = 0. 
For the Hubble parameter (1131) . the scale factor a{t) reads, 


> , fo(-t+t a ) 1+ ' 

a(t) = e !+“ 


(67) 


and therefore, a reads, 


/o(-t + t 3 ) 1 + 


/Q(-t + t s ) 1 + ° 


d(t) — e 1+ “ /o (—t + t s ) a — e 1+ “ fo{~t + t s ) 


\—1+a 


which can be recast as follows, 


( 68 ) 


a(t)=e /0< 1 +“’ fo(~t + t s ) 1+a (f 0 (-t + t s ) 1+a -a) . (69) 

Therefore, when a < 1, the parameter a never vanishes and is always singular, while for a > 1, it vanishes in two 
points, namely, 


t = t s 


/a\ 1/(a+1) 

V/07 


t = t s . 


(70) 


Notice that for a > 1 the physical system of the scalar field develops a Type IV singularity, as we demonstrated 
earlier. Let us elaborate on these cases, to see the physical consequences these bring along. For the case a > 1, if 
inflation ends at t = t s , the physical system has a Type IV singularity and a singular w e g at t = t s . with inflation 
ending at t = t s . if t s is considered to occur early enough. A more interesting case occurs if t s is considered to occur 
at late times, and inflation ends at a time t e , with t e , 


te — t s 


/ay/ (Q+1) 

\fo) 


(71) 


This can be true if fo a. In this case, inflation ends at i e , which might be quite earlier than t s , and at t s the 
Universe experiences a Type IV singularity. Interestingly enough, if inflation ends at the Type IV singularity solution 
t = t s given in Eq. GOD, which occurs if a > 1, then after t > t s , the effective equation of state describes quintessence 
acceleration. We have to mention, however, that, if we include quantum effects, the vacuum might acquire higher 
energy and consequently, inflation might occur again. So the analysis we performed is valid at a classical level, and 
in order to avoid such issues which originate from the quantum fluctuations of the scalar field, the hybrid inflation 
scenario has to be used fl5| . For the hybrid inflation consequences, the reader is referred to Ref. flbt . where this issue 
has been studied in detail. 
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Before closing this section, we shall briefly mention how the presence of perfect matter fluids affects the cosmological 
singularities of the physical system that consists of a scalar field with Hubble rate given in Eq. (fl3l) . In this case, the 
scalar field potential and u)(ip) become 0, 

w(</>) = - ~2 /'(<£) - (Wm + i)F Q e ~ 3{1+Wm)FW , 

V(4>) =4 ( 3 /^) 2 + /'(*)) + ^-^F 0 e-^ 1+w ^ F W , (72) 

where Fq an integration constant and F'(<j )) = From Eq. (1721) it easily follows that the exponential in both 

V((f>) and u>((/>), has a milder singularity structure, in comparison to the rest singular terms, for the Hubble parameter 
given by Eq. m Consequently, the same analysis we performed previously applies and in addition notice that 
the potential is not of a power law type, therefore it deviates from our primary motivation for our study, which is 
power-law potentials, so we omit this study for brevity. A detailed analysis on this issue was performed in Ref. jl5j | . 


II. UNIFICATION OF TYPE IV-SINGULAR INFLATION WITH DE SITTER DARK ENERGY EPOCH 

LEADING TO A TYPE IV FUTURE SINGULARITY 

As we already mentioned in the introduction, it is possible that the Universe may have passed from a Type IV 
singularity between the singular inflation era and the ACDM epoch. In addition, in the context of unified inflation- 
dark energy models given in Ref. [13], it is possible to have a future singularity of Type II, III or IV, because of the 
epoch of singular inflation, and in effect we may have a de Sitter-like or pure ACDM dark energy era. In view of 
these evolutionary aspects, in this section we shall provide an example, in which, we may have a unified description of 
singular inflation and quintessential dark energy epoch, leading to a Type IV singularity. The model we shall present 
is particularly fine-tuned to realize this scenario, nevertheless it is an example of a scalar-tensor theory leading to 
singular inflation and quintessential type of dark energy with a Type IV singularity at finite time. For a recent 
scenario in the spirit of ours, see [24| . 

We shall assume that the function defined in Eq. (fTUI) . is equal to, 

/(f) = Ci ( t + t/)“ + c 2 ( t + t s f + cit“ + g , (73) 

where g is an arbitrary constant. In addition, we impose the condition that the physical system of the scalar field 
described by Eq. (ESI), experiences a Type IV singularity at the time tf which is considered to occur at the end, or 
right after the inflationary era. 

In addition, it is also required that the Universe experiences an additional Type IV singularity at the time t s , which 
is considered to occur much more later than the present time, which we denote t p . So in effect, we have the following 
conditions, 


f s > t p > t f . 


(74) 


Thus the model (l73l) describes singular inflation plus a dark energy era, which we now study in detail. Since, the two 
singularities that occur are of Type IV, the parameters a and b must satisfy the following constraints, 


a > 1, /3 > 1. 


(75) 


and also we shall assume that a and b are integers. Using the reconstruction method of the previous section, the 
kinetic function w(</>) and the scalar potential V(<f> ) of the corresponding scalar-tensor theory are, 


>{ 4 >) = - 


2(-aci(t/-0) 1+a — pc 2 (t s — 4>) 1+/3 + ac\(j) 1+Q ) 


, rf ,, —aci(t f — </) 1+a - /3c 2 (t s - (j>) 1+/3 + aci0 1+ “ + 2 (g + a(t f - </)“ + c 2 (t s - <j>)P + c 1( / a ) 2 

V W>) =- ~2 -> ( 76 ) 


while the corresponding EoS defined in Eq. is equal to, 


Wett = 


1 2 (acit 1+01 — aci (—t + tf) 1+01 — (3c 2 (—t + t s ) 1+/3 ) 

3 (<7 T Cit a ci(—t + t/) a + c 2 (—t + f s )^) 


(77) 


Let us consider the case that a is equal to, 


a = 2m + 1 


( 78 ) 
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thus it is an odd integer. Near the Type IV singularity at £ ~ £/, the EoS (1771) behaves as follows, 

2 (acit _1+ “ - /3c 2 (-£ + £ s ) _1+/3 ) 

Wefi = — 1-o- , 

3 (g + c\t a + C 2 ( £ + £ s )^) 

since the following two terms vanish when £ ~ £/, 

ci (—£ + £/)“ ~ 0, aci(-£ + £/) _1+a ~ 0 , 


(79) 


(80) 


since a > 1. The expression given in Eq. CS) can be further simplified, by considering the fact that £ <C £ s , when 
£ ~ £/, and also that /3 > 1, in conjunction with the fact that — £ + t s — t s 1. So eventually, Eq. (1791) can be 
rewritten as follows, 


2 (aci£ 1+ “ — /3c 2 £ s 1+ 

Weff = -1--- -2 

3 \ g + cit a + c 2 tsj 

Notice the term appearing in the numerator of the fraction in Eq. m- This term is negative, because t s £ and 
(3 > 1. Therefore, the EoS of state is w e s > —1, for £ ~ tf, which describes quintessential acceleration. We now 
proceed to examine the behavior of the EoS for cosmic times corresponding to present time t p . So for £ ~ t p , since 
£>>£/, and owing to a being equal to an odd integer (1751) . the following two terms appearing in Eq. (1771) can be 
approximated as follows, 



(—£ + £/)“’ 1 ~ (—£)“ -1 = £“- x , (—£ + £/)“ ~ (-£)“ = -t a , 


where we used the fact that since a is odd integer (a > 1), a — 1 is even. Thus, the EoS can be written, 


w e f{ = -1 + 


2/3c 2 (—£ + £ s ) 1+/3 
3 (g + c 2 (-t + t s Y) 2 ’ 


(82) 


(83) 


which describes quintessential acceleration. By taking into account that we assumed £ s t p and also that £ is of the 
order £ ~ £ p , the equation above can be further simplified, to the following expression, 


W e ff = -1 + Ai, 


(84) 


where we have set, 


Ai 


2 pc 2 t~ 1+ P 
3 (g + c 2 ts^J 


= const. 


(85) 


By choosing g 1 and since £ s 1, in conjunction to the fact that /? > 1, the constant parameter is approximately 
zero, Ai ~ 0, and therefore the EoS becomes nearly w e g ~ —1, which is de Sitter acceleration. Notice that since we 
assumed £ ~ £ p , with t p corresponding to present time, this era is the dark energy, late time era. So it is possible to 
end up to a de Sitter dark energy era. 

Coming back to Eq. m, let us study how the EoS behaves near the future Type IV singularity that occurs for 
cosmic times of the order £ ~ £ s . Owing to the fact that £ ~ t s and also since £ £/, Eq. (1771) can be approximated 

as follows, 


w e ff = -1 


2 (aci£ 1+01 — aci(-t) 1+Q ) 
3(5 + ci £“ + ci (-£ + £/)“) 2 


since for £ ~ t s the following term vanishes, 


( 86 ) 


-/3c 2 (-£ + £ s ) 1+/3 ~ 0, 


(87) 


Owing to the fact that, for an odd a, we have (—£)“ -1 = £ a ^ 1 , the second term in Eq. (l86l) vanishes, so the EoS is 
approximately equal to w e s — —1. Consequently, near the Type IV singularity, the Universe accelerates in a nearly 
de Sitter way. In principle, the Universe could survive the future Type IV singularity, so we could also study the 
behavior after the future singularity, but we omit this, for brevity. We only highlight the possibility that the Universe, 
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after the future Type IV singularity can evolve to phantom acceleration, if /3 is equal to some even integer (which 
would imply that /3 — 1 is an odd integer). Indeed, the corresponding EoS, which for t > t s is equal to, 


w e ff = -1 + 


2(3c 2 (-t + t s ) 1+/3 
3(ff + c 2 {-t + t s )P) 2 ’ 


which is Weff < —1, since for t > t s and with /3 — 1 being an odd integer, the following holds true, 

2 /3c 2 (-t + t s )~ 1+l3 = —2/3c 2 (t - t s )~ 1+p . 


( 88 ) 


(89) 


Recall that we assumed a = 2 m + 1, and this was the most physically appealing case. The other case, a = 2 m yields 
similar results, but there is no possibility to have de Sitter dark energy, so we omit it for brevity, since the analysis 
is the same as in the case we just studied. 

Therefore, in this section we presented a scalar-tensor model, according to which, the Universe experiences the 
following: 

• A Type IV singular inflation 

• Possible quintessential or approximately de Sitter dark energy era 

• A future Type IV singularity after the dark energy era 

There are more possibilities to study, in the context of this model, for example if someone assumes that the future 
Type IV singularity occurs near the present time, that is t p cs t s . Also it would be interesting to find a pure modified 
gravity description with this type of behavior. We hope to address these issues in a future work. 


III. F(A)-GRAVITY ANALYSIS 

Having provided a quantitative analysis of the singularity structure in the context of singular inflation, with special 
emphasis in the Type IV singularity, we now study how the physical system under study can be generated by an 
F(R) gravity. Particularly, we shall be interested in which F(R) gravity generates the Type IV singularity, that is, 
how the F(R) gravity behaves near the Type IV singularity. For a similar approach to our reconstruction study, see 
Ref. [25]. Note that we can equivalently assume that the singularity occurs at the end or after inflation, or even at 
late times, since our approach is based on the existence of the Type IV singularity and also works as the cosmic time 
approaches the finite time where the singularity occurs. Therefore, having the Jordan frame F(R) gravity at hand, we 
can calculate the inflationary parameters to check the viability of the F(R) gravity, but we defer this to a future work. 
Before we start off, it is worth describing in brief the underlying theoretical framework. For reviews and important 
papers on the subject, the reader is referred to [H, H^]. The general F(R) gravity action, with matter fluids present, 
is given by, 


S =^2 f d 4 Xy/^F(R) + S m (g pv , T m ), (90) 

where n 2 = 87 tG and in addition, S m describes the matter fluids action. We shall adopt the metric formalism approach 
[l(|, in the context of which, by varying the action with respect to the metric g pv , we obtain the following equations 
of motion, 


Ku - \r. 9 ,, v = (t ijm + 1 f F(R) 2 RF,(R) g^ + V m V„F'(R) - g^DF'iR)^ ) . (91) 

with the prime denoting in this case differentiation with respect to the Ricci scalar R and moreover, T pu stands for the 
energy momentum tensor corresponding to the matter fields. In the F(R) modified case, the energy momentum tensor 
receives an extra contribution originating from the F(R) gravitational sector, with this addition being the feature 
that renders F(R) gravity a modified version of ordinary Einstein-Hilbert gravity. This additional contribution to the 
energy momentum tensor is equal to, 


rriefi 


1 l'F{R)-RF'(R) 
k V 2 


+ V^Vj ,F'(R) — g^OF’ (R) 


K 


2 


(92) 
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Finally, in all the following sections, we shall assume a flat FRW metric of the form given in Eq. ©, in which case 
the Ricci scalar is given by, 

R = 6(2 H 2 + H ), (93) 

Before providing a full quantitative description of the F{R) gravity near the singularity, we shall give a qualitative 
description of how the F(R) gravity behaves near the singularity. We assume that the F(R) gravity is of the form, 

F(R) ~ F 0 + F x R e . (94) 


where Fq and F\ are constants, with Fq an arbitrary constant, allowed to be even zero, but F\ is constrained to be 
non-zero, that is, F\ ^ 0. Then by using the trace part of the F(R) equation, 


3D/'(R) =R + 2f(R)-Rf'(R). 


(95) 


we obtain the following relation, 


3FiDf? e_1 = 


R 


when e < 0 or e = 2 


(2 — e) FiR e when e > 1 or e ^ 2 
Then in the FRW background (|5Jl. the Hubble rate H behaves as follows, 

TT _ k 0 

±1 a « 

(t s - t) P 

where ho and /3 are constants. Consequently, the scalar curvature R = 6 H + 12 H 2 behaves as follows, 

12 hi 


(96) 


(97) 


(ts-t) 


w 


when /3 > 1 
R ~ { when P = 1 

when /? < 1 


(98) 


According to Eqs. (l97l) and (j98|) . it easily follows that the /3 > 1 case corresponds to the Type I (Big Rip) singularity 
in [id, [l2f, and the 1 > f3 > 0 case corresponds to Type III singularity. Moreover the 0 > (3 > — 1 case corresponds 
to Type II while when /? < — 1 (but /3 ^ integer) we have a Type IV singularity. 

We now proceed to the detailed description of the F(R) gravity near the Type IV singularity of Eq. (fl3l) . by using 
the reconstruction techniques that were firstly developed in |27H29j . What we are mainly interested in, is to find 
which F(R) gravity exactly generates the Hubble rate (fl3l) . for /i(t) = 0 and f 2 {t) = fo, as t tends to t s . Notice that 
as t —> t s , then t — t s —> 0. The action of a general F(R) gravity, in the absence of matter fluids, is equal to, 


S = ^f d A xV^F(R). (99) 

By varying the above action with respect to the metric, we obtain the first FRW metric which is equal to, 

- 18 (AH{t) 2 H(t) + H(t)H(t)^ F"(R) + 3 ( H 2 (t ) + ij(t)) F\R) - = 0 . (100) 

Using an auxiliary scalar field <f>, the F(R) gravity action given in Eq. (1991) can be written in the following way, 


S = j d a xV^9(PWR + QW) ■ 


( 101 ) 


The main goal of the reconstruction technique is to find the exact form of the functions P((f>) and Q((f >), which depend 
explicitly on the auxiliary scalar field. Particularly, if these functions are expressed as implicit functions of the Ricci 
scalar, the F(R) gravity will follow easily. In order to find these, we shall use the fact that the auxiliary scalar field 4> 
is a dynamical (meaning that it depends on time) degree of freedom, auxiliary because it has not kinetic term. Then, 
by varying the auxiliary field modified action of Eq. (11011) with respect to this auxiliary degree of freedom, we get the 
following algebraic equation, 


P'{4>)R + Q'{<j>) = 0 


(102) 
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with the prime this time denoting differentiation with respect to <fi. If this equation can be solved explicitly with 
respect to the auxiliary field </>, will yield the function <f){R). Correspondingly, the F(R) gravity can easily be obtained 
by making use of <j>(R) and substituting this to the auxiliary F(R) action of Eq. 11011) . and we finally get, 

F(4>(R))=P((KR))R + Q{<f>(R)). (103) 


Hence, what we need to do is to find the functions P(<f>) and Q((f >), given the Hubble rate, 

H(t) = f 0 (-t + t s ) a . 


(104) 


In order to find the functions P(4>) and Q(4>), we vary the action (1 101 D with respect to the metric tensor, and by using 
a flat FRW metric, we obtain the following differential equation, 


0 = — 6 H 2 P{(j)(t)) - Q{4>(t)) - 6 H 


d pm)) 


= 0 , 


0 = (AH + 6P 2 ) P(0(t)) + Q(<£(£)) + 2 


d t 

d 2 pm)) , d pm)) 


d t 2 


At 


= 0. 


By eliminating Q(<j>(t)) from Eq. (11051) we obtain, 

0 d 2 pm)) 


At 2 


- 2 H(t )+ 4 HP{<j>{t)) = 0 . 


(105) 


(106) 


Therefore, by solving Eq. (11061) with respect to P(£), we obtain the explicit form of P(t), given the Hubble rate. 
Accordingly, we easily find the function Q(t), by using the first equation of Eq. (11051) . Notice the interplay of 
the variables <j> and t, as they appear in the previous equations. As is explicitly proven in Ref. [27j . owing to the 
mathematical equivalence of the F(R) action (l99l) and the auxiliary F(P) action 11011) . the scalar field is identical to 
the cosmic time t , that is </> = t. 

Using this technique, substituting the Hubble rate 11041) in Eq. 11061) . we get the following differential equation, 


- 2/o(—£ + t s ) a ^- - 4/o(—£ + t s )~ 1+a aP(t) = 0. 
By setting x = — t + t s , Eq. (11071) can be written in the following way, 


d 2 P(%) 

Ax 2 


dP(x) _ 2/ -i+« aP(;r ) = o . 
ax 


(107) 


(108) 


This differential equation can be easily solved analytically by substituting z = x a+1 , in which case the differential 
equation becomes, 


(a + 1 ) 2 z —+ (a + l)(foz + cz )—^ — 2foaP(z ) — 0 , 


with general solution, 


P(z) 


ln(z) — z/q 

= e l+« 


CiU 


-1 - 3a 1 zf 0 \ 

1+a ’ 1+a 1+aJ 


ln(z)-z / 0 

+ e !+« 


C 2 L 


m 

n 


zfo \ 

1 + a) ’ 


where the parameters n and m are equal to, 


n = 


-1 - 3a 
1 + a 


m = 


1 

1 + a ’ 


(109) 


( 110 ) 


( 111 ) 


also C'i,C 2 are arbitrary constants and finally U(a,b,z) and L™(z) are the confluent Hypergeometric function and 
the generalized Laguerre polynomial respectively. By substituting again z = x a+1 we get, 


P(x) 


= e r+s 


CxU 


-1-3a 1 x a+1 fp \ 

1 + a ’ + 1 + a’ 1 + a ) 


ln( I ° + 1 )-x°+ 1 .f, ) 

+ e 


C 2 L 


m 

n 


(t^) 


( 112 ) 
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Having at hand P(x ), we can easily compute Q(x), which reads, 


Q(x) = - -^—6e-^T^f oX ~ 1+a (x 1+a ) ^ ((1 + f 0 x 1+a ) 


1 + a 

x (1 + a) ( CiHl 1 


1 x 1+a fo \ , r rmi 

1 + a ^ ' 1 + a’ 1 + a > ° 2 ni 


x ( (1 + a) C\U ( 4±i^,l 


1 x 1+a f 0 


9 + 4 a 

+ (l + 3a)Cif/( +—1,2 


1 + a’ 1 + a’ 1 + a 
1 x 1+a fo 


1 + a’ 1 + a’ 1 + a 

x 1+a fo 
1 H- ot. 


+ (1 + a) C' 2 + K 2 2 


fo 


X 1+a fo 
1 -|- Gi 


— X 


1+a 


(113) 


(114) 


with rii, mi being equal to, 


m = -■ 


1 + 3a 
1 + a ’ 


mi 


1 

1 + a ’ 


while ri 2 and m 2 are equal to, 


n 2 = —■ 


2 +4a 
1 + a 


m 2 = 1 + 


1 + a 


(115) 


(116) 


Substituting P(x) and Q(x) from Eqs. (11121) and (11131) in Eq. (11021) and solving with respect to x, one can in principle 
derive the final form of the function x(R). This task, however, is rather formidable, owing to the complex structure of 
the functions P{x) and Q(x). In order to obtain an approximate form of the resulting F(R) gravity, we shall exploit 
the fact that we are interested to know the F(R) gravity near the Type IV singularity, so this means that x —> 0. We 
therefore expand P'(x) and Q'(x) in small powers of x, and by keeping dominant terms (recall that a > 1 in order a 
Type IV singularity occurs), the resulting expansion of P(x) reads, 

P ^^T+°( a; )’ ( 117 ) 


with the parameter A being equal to, 


c > r (+) (ife) -1 * 5 

(l + 3a)C 1 r(l + T ^) 

(l + a)r(^) 

( i + a)r (i + it^) 


Moreover, the function Q(x) for small values of the argument x can be approximated as, 

Q( x ) — xa+1 + c + o{x ), 


(118) 


(119) 


with the parameters B and C being equal to, 


6 / 0 (l + 3a)Ci | 

((2 + a)T| 

( 3+5+ 

I r | 

( 1 + T^. 

) -2(1 + 

2a)T | 

( 2+4a ^ 

^ !+“ J 

l r i 


0 

(i+ a ) 

,-T+ 

(l + a) 2 r( 

2+4a ^ 

. !+“ J 

I r | 

f 3+5 a^ 

l 1+a j 

) 


/'(l+a)(2+a)r( T k) Q (l+3a) 2 r(l + T +j) 4 (l+5a+6a 2 )r(l + T i J ) , 4(H-2a) 2 (l+3a)r(2+ T + J ) 

J ° 1 ^ ^-m) r(3jg) r(^) + (2+a)r(^) 

(1 + a) 3+ “ 


( 120 ) 


Substituting Eqs. (11171) and (11191) in Eq. (11021) . we finally get, 


x + - 


AR + B 


( 121 ) 
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Finally, by substituting (11211) in Eq. (11031) . we obtain the final form of the F(R) gravity that generates the Type IV 
singularity, which is, 


A 2 RA R 2 

F{R)^-^R 2 -2-^R- — +C. (122) 

Notice that this F(R) gravity can be brought to a form of Einstein-Hilbert gravity plus curvature corrections, if the 
free parameter C\ is chosen appropriately. Indeed, if, 



( 123 ) 


which can be true if C\ is equal to, 


C(l + a) 2 r( 

f 2 + 4 a 

, i+“ J 

1 T | 

( 3 + 5 cx' 

l !+“ , 

) 

2 A 6/ 0 (1 + 3 a) | 

((2 + a)T| 

( 3 + 5 a) 

l !+“ ) 

l r i 

( 1 + ii 7 y 

) - 2(1 + 

2 a)T( 

f 2 + 4 a^ 

, i +«) 

|T| 

+ +))(*) 



then, the F(R) gravity is, 


( 124 ) 


A 2 R 2 

F(R)~R-—R 2 - — +C. 


( 125 ) 


Thus, the resulting F(R) gravity near the Type IV singularity can be brought to an Einstein-Hilbert form, plus 
curvature corrections and cosmological constant. Notice that the final form of the F(R) gravity near the singularity 
is of the form (1941) . that is, of polynomial type. 

Actually, the form of F(R) gravity is a nearly R 2 gravity plus cosmological constant, when 4r- < 0. In order to see 
this, we change the notation of the coefficients in Eq. (I125D as follows, 


C B 2 

C » = -4X' A = “ +C ' 


By doing so, the Jordan frame F{R) gravity takes the following form, 


( 126 ) 


F(R) = R+§ r + A. 
4Co 


( 127 ) 


This Jordan frame theory is equivalent to a variant of the Starobinsky model 21] in the Einstein frame. The pure 
Jordan frame F(R) gravity is of the form, 


S = ^~>j ^V^F(R), ( 128 ) 

Starting from action (11281) . and by introducing the auxiliary field A , the action (11281) can be cast in the following 
form, 


s= ^S + f ^ ■ ( i29 ) 

By varying (I129D with respect to the auxiliary field A, we easily obtain the solution A = R, and consequently this 
verifies the mathematical equivalence of the actions (11281) and (11291) . We can be transferred to the Einstein frame by 
making the canonical transformation, 


V = ( 130 ) 

where tp is the canonical Einstein frame scalar field (or inflaton field). By conformally transforming the Jordan frame 
metric, 


9^ = 


6 ^ 9 ^ ! 


( 131 ) 
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with the “hat” denoting the Jordan frame metric, we get the Einstein frame canonical scalar field action, 

2 




ni ,M d Ad A _— ( A _ 

g d„Ad„A 2fc2 F ^ A y 


y 2 k 2 2 Vf '(^)7 

= J d i x y /=g - ^g^dptfidvip - V . 

The scalar potential V{tp), which is a function of the canonical scalar field ip is equal to 
A F(A) 1 


(132) 


v(v) = 


\/2k 2 /3¥>^ _ g 2- v /2K 2 /3 <f F F ^ e ~y/2K 2 /3<f'j ^ _ ( 133 ) 


F'{A) F'{A) 2 2 K 2 

Upon using the form of the Jordan frame F(R) gravity (11271) . the scalar potential of Eq. (11331) becomes [H], 

V(tp) ~ C 0 + C 2 e~ 2 ^ KV + Cie-V^I^ , (134) 

with the coefficients C\ and C 2 standing for, 

Ci =-2C 0 , C 2 = C 0 -A. (135) 

The scalar potential (I134D describes a nearly Starobinsky Einstein frame canonical scalar theory. If the coefficient C 2 
is equal to Co, then the potential (I134D becomes exactly the Einstein frame Starobinsky model, namely, 


U(^) = C 0 (l-e\/i ^) 2 . 


(136) 


IV. IMPACT OF TYPE IV AND FINITE TIME SINGULARITIES ON SCALAR PERTURBATIVE 

MODES 

In this section we shall briefly discuss the impact of finite time singularities on the power spectrum of perturbed 
modes, with special emphasis to the Type IV singularities. A thorough analysis on this issue however was presented 
in (30j| . with regards to the classical stability of sudden singularities. In addition, in f30f the Big Rip case was also 
extensively discussed. Notice that the Type IV singularity falls in the category of some sudden singularities and before 
we proceed to a qualitative analysis of the model we studied in this paper, it is worth outlining the most important 
results of [30]. For details on this important issue, the reader is referred to [30( and related references therein. As was 
shown in |30[, if a(t) ~ (t s — t) x , the energy density and pressure are equal to p ~ (t s — t) 2X ~ 2 and p ~ (t s — t) x ~ 2 . 
Then if A > 1 and if A is a non-integer number, all the scalar metric perturbations are bounded. Let us make some 
correspondence of the finite time singularities with sudden singularities. The case 1 < A < 2 describes the Type 
II singularity, so the results of [30] indicate that singularities like Type II or milder can provide us with bounded 
scalar metric perturbations. Actually, the Type IV singularity is milder than the Type II, so it is expected that the 
scalar metric perturbations are bounded in the Type IV case too. In the rest of this section we shall perform a brief 
qualitative analysis for the case at hand, with special emphasis to the long-wavelength modes and for a Type IV 
singularity. 


A. Long-wavelength solutions during and after inflation: The Type IV case 

As is very well known L, f3ll . 13^ , the quantum fields which are present during the primordial acceleration of the 
Universe, experience quantum mechanical fluctuations. These primordial fluctuations could serve as a very crucial 
observational feature and a test of any potentially viable cosmological theory. This is owing to the fact that the 
primordial perturbations could have a direct impact on the present epoch Universe’s large scale structure. In principle, 
the kinetic energy of the scalar fields is responsible for the generation of a density distribution at a given spacelikc 
hypersurface, during the primordial acceleration. The perturbations of this density distribution are responsible, to a 
great extent, for the large scale structure in the Universe, materialized by gravitational instabilities. 

Consider a homogeneous, isotropic and spatially flat Universe. The quantity that the determines the scale that 
within it causal physical processes can be materialized, is the sound Hubble radius c s H _1 , where H is the Hubble 
rate. The parameter c s stands for the sound speed, which for a canonical scalar field c s = 1 (this corresponds to 
the case we study in this paper). The sound Hubble radius plays a crucial role towards the understanding of the 
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way that the spectrum of primordial perturbations is actually generated. Another physically important quantity is 
the fraction k/aH 1 with a the scale factor and k the wavenumber. As the primordial acceleration occurs, the scale 
factor increases rapidly, and consequently the comoving Hubble radius ( Ha ) _1 , decreases. During the same epoch, 
the perturbation’s physical length , satisfies the relation a/k -C c s Jf _1 . The latter relation is very important since this 
actually determines the way that the scalar modes are evolving. We shall quantify this evolution later in this section, 
but let us qualitatively analyze what is the meaning of this. Practically, the relation a/k <C c s 7? _1 indicates that the 
physical scales grow more quickly in comparison to the Hubble sound radius and in effect, the scalar modes we study, 
begin their evolution within the horizon c s i? _1 . In the process of the accelerating expansion, the perturbation modes 
(scalar field quantum fluctuations), become larger more quickly than the sound Hubble radius. Consequently at some 
point it passes out of the horizon, so it freezes out. This freezing has a direct impact on the scale dependence of the 
scalar perturbations, at the time this freezing occurs, which is when the Hubble sound horizon is crossed. Practically, 
the modes at that point enter the long-wavelength limit aH 3> c s k. In the context of this approximation, what plays 
a crucial role for the determination of the scale dependence, is the equation of state of the matter fluids present at 
moment of the horizon crossing, which can occur for example during the matter or during the post recombination 
era. In the following we are not interested in determining the evolution of the mode per se, for a specific matter fluid 
present, but our aim is to study the effect of a Type IV singularity on the evolution of a long-wavelength frozen mode, 
with emphasis to the question whether the modes are bounded or not. 

Assume that the singularity appears at the end of inflation, so t s indicates the time that inflation ends. The scale 
factor corresponding to the Hubble rate of Eq. m, with h(t) = 0 and f 2 (t) = /o, namely H(t) = f 0 ( t s - t)“, is 
equal to, 


a(t) 


/n(—*+ t s) 1+l 

= e i+“ 


(137) 


What is the main focus in this section is to study the behavior of the scalar perturbations during and after the 
inflationary era and qualitatively check whether the existence of a Type IV singularity at the end of the inflationary era 
affects the spectrum of primordial scalar perturbations. The perturbed metric that will yield the scalar perturbations, 
is of the following form, 


ds 2 = (1 + 2<p) d t 2 — a 2 (l — 2ip) da; 2 , (138) 

i 

where <p and ip quantify the metric perturbation. Notice that in order to obtain the scalar field perturbation, we also 
perturb the canonical scalar field ip, in the following way, tp(t,x)' = (p + d<p(t,x). Note that x denotes all the space 
coordinates and (p denotes a homogeneous scalar field. Since at first order in perturbation theory, <p = ip, in order 
to obtain the metric perturbations only, we are left practically with two variables, the variable <p and the inflaton 
variation 5<p , which are related through the Einstein equations. It can be shown [3lj that the equation which describes 
the evolution of one scalar mode <p k , with wavenumber k, takes the following form, 

<pk +(h- 2^j 4> k - (2 H% + 87tG+ - <p k = 0, (139) 

with a(t) being the scale factor. We shall be interested in the long-wavelength solutions of this equation, which 
corresponds to the case aH c s k. Since the speed of sound for a canonical scalar field is equal to one, the limiting 
case at hand corresponds to wavelengths with aH k. It is possible to show that during any stage, the long- 
wavelength (fc —> 0) analytic solution of Eq. (11391) takes the following form, 

cp k = Gr(fc) (l - ^ J* a{t)dtj - > ( 14 °) 


with Ci(k) the coefficient of the non-decaying modes, while C 2 (k), the coefficient of the decaying modes. Using the 
scale factor of Eq. (I137D , we easily obtain the following form of the long-wavelength mode (p k , 
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It is obvious that the existence of a Type IV singularity does not cause any inconsistencies to the long-wavelength 
modes, with inconsistencies meaning any type of singular or abnormal behavior. Therefore we may conclude that 
the Type IV singularities are harmless. A qualitatively similar conclusion was derived in [§] (see also (30j| ). with the 
corresponding transition to the standard cosmological model not affecting the power spectrum. 


V. BRIEF DISCUSSION OF FINITE TIME SINGULARITIES AND THE GRACEFUL EXIT OF 

INFLATION 

One of the most serious problems that the first models of inflation [33!] were confronted with, was the exit from 
the inflationary era, a problem known as the graceful exit. Initially, in the context of the false vacuum theories 
of inflation, it was not possible to consistently describe a mechanism for the smooth transition of an inflationary 
Universe to the classical FRW expansion. This problem was solved in the new inflation theories [U, [3S| . where the 
slow-roll approximation actually solved the problem of graceful exit. Indeed, if the inflaton scalar slowly rolls down the 
potential, towards the minimum of the potential, then inflation ends where the slow-roll approximation breaks down, 
and the field oscillates around the minimum of it’s potential. In the case we studied in this paper, the canonical scalar 
field potential of Eq. (1281) belongs to a class of potentials were slow-roll applies. Therefore, the graceful exit issues 
in our case is dealt within the context of power-law inflation models, which are in good agreement with observations 

G3- 

However, the exit of inflation in our case becomes a serious problem when the slow-roll approximation breaks down. 
As we demonstrated in a previous section, this is true when the Type IV singularity occurs during the inflationary 
era, since the slow-roll parameters become divergent. Also the e-folding number becomes divergent at the time the 
singularity occurs. Therefore, the problematic case is only when the singularity occurs during the inflationary era. 
If however the singularity occurs at the end or much later than the inflationary era, the slow-roll approximation is 
not violated and therefore the graceful exit can be dealt in the context of the power law models of inflation. If the 
singularity occurs much more later than the ending of inflation (even during the dark energy era), then slow-roll 
breaks when the inflationary indices become of order one, and the exit from inflation goes in the standard way, for 
example, the resulting de Sitter Universe becomes unstable and the exit is guaranteed. 

Before we close this brief discussion, it is worth to note another quite physical appealing but exotic new possibility 
that a singularity during inflation might generate. As we demonstrated in the text, in the power-law models of inflation 
of the form (1281) . when the singularity occurs during the inflationary era, the slow-roll approximation is violated in a 
violent way, since the slow-roll parameters become divergent and also the e-folding number it self becomes divergent. 
If however the violation of the slow-roll approximation is not divergent, this could be viewed as another mechanism 
for a graceful exit from inflation. So the appearance of a singularity during inflation could generate instability in 
the inflationary spacetime, and this could indicate a new mechanism for a consistent exit from inflation era. This is 
another view to the Type IV singularities and work is in progress towards this direction. 


Conclusions 

We have thoroughly examined the existence of classical finite time singularities in the context of scalar-tensor 
theories. Using a very generic behavior for the Hubble rate, we provided a quantitative description of all the finite 
time singularities of the physical system under study, with emphasis to the Type IV singularity. Particularly, with 
regards to the Type IV singularity, we investigated when this can occur and we also showed what is the impact of this 
singularity on the slow-roll inflation parameters. As we explicitly demonstrated, the presence of this type of singularity 
may have dramatic effects on the slow-roll parameters, rendering them singular, in the case that the singularity occurs 
during the inflationary era. After exemplifying our claims with illustrative examples, we studied the behavior of the 
effective equation of state parameter and in which cases, its value implies acceleration. In addition, we thoroughly 
examined all the cases at which the Universe might experience a finite time singularity, with the singularity being 
during, at the end, after inflation. In addition, we studied the case at which the singularity occurs much more later 
than inflation ends. The issue of when inflation classically ends was examined too. As we showed, in the case that a 
Type IV singularity is present, the results are of some interest, especially in the particular case that the time t s , at 
which the singularity occurs, is much more later in comparison to the time t e , at which inflation ends. 

In the context of the above scenarios, the Universe experiences a Type IV singularity at finite time t s . Interestingly 
enough, the Type IV singularity is not a crushing singularity, such as the Big Rip singularity is for example. This 
means that at that point, extensibility and geodesic incompleteness do not necessarily occur. With regards to this, 
the nature of the Type IV singularity needs yet to be understood at a fundamental level. Usually in nature, non 
extensible singularities are protected, as in the case of spacetime singularities in black holes which are dressed with 
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the horizons of black holes. So having at large scales undressed, naked, crushing singularities, is a rather inconvenient 
feature of classical cosmological theories. Perhaps the crushing singularities occur as a result of our classical approach 
and maybe it is the signal that a quantum gravitational description is needed to amend these issues in a concrete 
way. In fact, in the context of loop quantum gravity, initial crushing singularities do not occur, see for example 
[36f |. However, the Type IV singularity is not a severe spacetime singularity, so it is of critical importance to find a 
consistent explanation for the occurrence of unbounded values of the higher derivatives of the Hubble rate, an also 
critically examine if this feature can be considered harmful for the theory itself. A first step towards this direction 
was attempted in this paper, where we showed that if the Type IV singularity occurs during inflation, this can be 
catastrophic for the slow-roll parameters and the corresponding observational indices. But if this Type IV singularity 
occurs at the end, or right after the end of inflation, compatibility with the observational data can be achieved to 
some extent. 

Another feature that has to be addressed, in the context of finite time singularities, is the issue of including quantum 
corrections, in the spirit of Ref. i9|. For the case of the Type IV singularity, quantum effects would be negligible, owing 
to the fact that the curvature is not growing significantly there. However, the quantum effect of conformal matter 
around the Type I, II and Type III singularities should be taken into account. As was shown in |sj, the singularity 
structure gets milder in the presence of quantum corrections, so the same is expected in the present case too. 

Finally a comment is in order. Having the possibility of a singular inflation, makes in some way the need for 
more accurate observational data more compelling. This is due to the fact that in this way, the various choices of 
a consistent theoretical description of inflation will eventually be narrowed down. With the present paper, we have 
excluded scalar-tensor models which cause a Type IV singularity during inflation, since these would be in conflict 
with current observational data, and also we have set the stage for new investigations towards this line of research. 
Therefore, in the absence of an ultimate theory of quantum gravity that will explain all the singularity issues, even a 
better understanding of non-crushing singularities, like the Type IV one, is of critical importance. 
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